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We consider the vacuum energy for a configuration of a sphere in front of a plane, both obeying 
conductor boundary condition, at small separation. For the separation becoming small we derive 
the first next-to- leading order of the asymptotic expansion in the separation-to-radius ratio e. This 
correction is of order e. In opposite to the scalar cases it contains also contributions proportional to 
logarithms in first and second order, elne and e(lne) 2 . We compare this result with the available 
findings of numerical and experimental approaches. 
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I. INTRODUCTION 



The Proximity Force Approximation (PFA) is the most important approximation for calculation of forces between 
curved surfaces at small separation. It originates from a work on adhesion back in 1934 based on the simple idea to 
integrate the local force density known from parallel surfaces [l[ . Given a sufficiently fast decrease of the force with 
separation this method works independently from the kind of force. In this respect it is universal. However, this 
method does not allow to beyond nor does it give any information on its precision. Attempts to repeat the original 
idea 'plane based' or 'sphere based', turned out to give misleading results. It was only with the new method of 
calculating the Casimir force, the T-matrix or TGTG-formula approach that the door opened to go beyond PFA. In 
these approaches, the interaction energy is represented by a infinite dimensional determinant. At large and medium 
I , separation this matrix can be truncated to become finite (even low) dimensional and numerical evaluation is possible. 
■ At small separation this does not work and below 

: , 
p; 

(d-distance, i?-radius of curvature, both at closest separation) the numerical effort is unmanageable. 
I ■ It must be mentioned that only e < 0.1 is the experimentally interesting region for forces between macroscopic 
^\ I bodies. This is because of the van der Waals and the Casimir forces as being quantum effects are generically mi- 
. croscopically small and become measurable only when multiplied by a microscopically large interaction area. Since 
the forces decrease proportional to d~ 4 at large separation, only the combination of small separation together with 
large radius of curvature allows for measurements with appreciable precision. A typical value used in experiments on 
^ ■ precision measurements of the Casimir force [2) is e ~ 10 -3 . 

The interest in corrections beyond PFA is triggered from both, theoretical and experimental sides. The first follows 
from the challenge to improve a situation which lasted more than 60 years, the second from the high precision of 
the contemporary force measurements and from the accuracy one would like to achieve for their comparison with the 
theoretical predictions. It must be mentioned that this has implications far beyond the atomic or solid state physics 
as mean to obtain stronger constraints on new physics [Fifth Forces), see for example |3[. 

For scalar fields, analytical corrections beyond PFA were obtained as an asymptotic expansion 

E 

l + ae + ... (2) 



-EpFA 



of the energy for e — > and simple numbers were obtained for the coefficient a. In Q this was done for the geometry 
of a cylinder in front of a plane for a scalar field obeying Dirichlet or Neumann boundary conditions. This includes 
the electromagnetic field at once since its polarizations separate in cylindrical geometry. The corrections for a sphere 
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in front of a plane were obtained in [5] for a scalar field, again for both, Dirichlet and Neumann boundary conditions. 
However, in this case the corrections for the electromagnetic field do not follow because of the non-separation of the 
polarizations. 

It is the aim of the present paper to fill this gap and to obtain the first beyond-PFA corrections for the electromag- 
netic field obeying conductor boundary conditions in the geometry of a sphere in front of a plane. Surprisingly, the 
asymptotic expansion in this case turns out to contain logarithms, i.e., it has the form 



The coefficients are calculated below and take the values a — —5.2, f3 — —0.0044 and 7 = 8.5 10~ 6 . The logarithmic 
terms come in from contributions which are specific for the vector case. 

In an experimental effort was undertaken to measure the corrections beyond PFA by using several spheres whose 
radii varied from 10 to 150 fim at separations d — 200 . . . 800 nm. The expansion was assumed to have the form of 
(|2"|) and the coefficient a was found to be zero within the experimental precision. Also numerical efforts are reported 
(for a cylinder in front of a plane in 0] and for a sphere in *8j, ||) by pushing the truncation in the T- matrix approach 
to higher orders and extrapolating towards the known value at zero separation. The results show agreement with 
the analytical results for Dirichlet boundary conditions but not for Neumann boundary conditions; details will be 
discussed in the last section. For a scalar field obeying Dirichlet boundary conditions results were obtained in [l(| 
using the independent method of world line approach. Like the extrapolation these confirm the analytical results. 

It is a second aim of the present paper to discuss in detail the analytical corrections beyond PFA for all combinations 
of boundary conditions for the scalar field. For instance, it will become evident that and why the corrections for a 
sphere in front of a plane are the same with Dirichlet boundary conditions on the sphere, but Dirichlet or Neumann 
boundary conditions on the plane. This case is interesting since the numerical results are different for the two cases. 

The paper is organized as follows. In the next section we give a representation of the vacuum energy in the sphere- 
plane geometry for all boundary conditions with special emphasis on the translation formulas used. In the third 
section we re-derive the asymptotic expansion for the scalar field and in the fourth section we derive the expansion 
for the electromagnetic field. In the last section we discuss the results. 
Throughout the paper we use units with h = c = 1 

II. REPRESENTATION OF THE VACUUM ENERGY IN SPHERE-PLANE GEOMETRY 

The representation of the vacuum energy in sphere-plane geometry was first derived in (llj | for the scalar case within 
the multiple scattering approach. Subsequently there appeared numerous variations of the derivation; we use that in 
Chap. 10. Thereby we highlight one essential step - the use of the translation formulas - having in mind their 
importance for understanding the logarithmic contributions in the electromagnetic case. 

The general structure of the T-matrix representation of the vacuum interaction is 




(3) 




(4) 



The geometry is shown in Fig.l. 




Fig.l: The configuration of a sphere in front of a plane. 
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In the frequency was rotated towards the imaginary axis, to — > i£. The symbol G denotes the Greens function 
(or, strictly speaking, the corresponding operator) describing the propagation from the sphere to the mirror and back 
and T is the T-matrix operator for the scattering on the sphere. Using the basis 

Uk,i m ( r ) = ji(kr) Yi m (Cl r ), (5) 

where the Yi m (Q r ) are the spherical harmonics and j(l(r) are the spherical Bessel functions, the Greens function can 
be written as 

Ira 

The limits of the summations are I > and |m| < I. 

In ([6]) both spatial arguments, r and r', appear to be in one and the same coordinate system with spherical 
coordinates (r, f2 r ). However, in the considered geometry, we would like to expand the T-matrix operator in its own 
coordinate system centered in (0, 0, d + R). If this system is taken for r, we need for r' a translation from this one to 
the mirror and back. This can be achieved by the translation formula 

Uk,i m (r + ae z ) = Ai mt i> m > (a) u k ,i> m ' (r) (7) 

I'm 1 

with a — 2L and e z is the unit vector along the z-axis. In (0, Ai m ^ m i(a) are the translation coefficients. These 
involve the Clebsch-Gordan coefficients, for details and an explicate expression see, for example, Eq.(10.125) in [l2T |. 
Applying these formulas, after some transformations, the energy (01 takes the form 

E=- [ ^Trln(l-iV), (8) 
2 loo 2^ 1 ' W 

where the trace is the orbital momentum sum, 

oo oo 
m=0 l=\rn\ 

In ([5]), N is an infinite dimensional matrix in the orbital momentum index, Nu>. Using the known relation Trln = 
lndet the energy can represented as a determinant. However, we do not use this in the following. 

The entries of the matrix iV are, for Dirichlet boundary conditions on the sphere and with the notation iV ; D ; , for 

— l+l' 



The function 



N %' = \l4£L £ Kl»+iM2ZL)H l u,dY(ZR). (10) 
^ l"=\l-V\ 

1 ^ ' K l+1/2 (£R) V ' 

is up to a factor the T-matrix for the scattering of a scalar field on a hard sphere in orbital momentum representation. 
The corresponding expression Nf 1 ,, for Neumann boundary conditions on the sphere can be obtained with 

dfm = (W£fl)/v5E)' (12) 

(K l+1/2 m/vm 

in place of dP(£R) in dTUJ) . In the above formulas I v and K v are the modified Bessel functions. We note that in 
opposite to eqn. (10.140) in the function I1+1/2HR) carries the index / in place of /'; a substitution which is 
allowed under the trace in (HI). In (fTTJl) we used the notation 



H l u, = V(2J + 1)(2J' + 1)(2Z" + 1) 
I V l"\ ( I V I" 



J \ m — ui 



(13) 
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which involves the Clebsch-Gordan coefficients coming in from the translation formula. We use the notation of the 
3j-symbols. 

Eq. {8]) represents the energy with Dirichlet boundary conditions on the plane. The case of Neumann boundary 
conditions on the plane is obtained by changing the signs in front of N. We unite all four combinations of boundary 
condition in 

£;XY = ^/li Trln ( 1 - ( - irivY )- (14) 

The first index denotes the boundary conditions on the plane, X — D with x — +1 for Dirichlet and X — N with 
x = — 1 for Neumann conditions. The second index denotes the boundary conditions on the sphere with Y = D for 
Dirichlet and Y = N for Neumann conditions. In the following we write all formulas for Dirichlet boundary conditions 
on both and discuss the other cases at the end of the next section. 

The logarithm in is taken of a matrix. In the following we will use the expansion of this logarithm such that 
the formula for the energy takes the form 

i r A -i 

E = -J ^£^-r (is) 



2 J m 2tt ^ s ■ 

J ~°° s=0 
oooo/soo\/s \ 

x£ e n £ n*wj 

m=0 l = \m\ \j=ll f = \m\J \i=0 I 

with the formal setting Iq = l s+ i = I. It should be mentioned that the Nij> are diagonal in the azimutal index m such 
that the sum over m appears only once. 

The expansion ()15|) corresponds to a perturbation series with respect to the T-matrix of the scattering on the 
sphere. This scries appears to converge. The convergence can be easily seen for large separations where it corresponds 
to the multipole expansion and also for short separations as we will see below. It is known that the vacuum energy 
for a transparent sphere, or for a background potential, where it makes sense to introduce a coupling constant, has 
the same structure as Eq.®. The last equation, in this way, appears as a perturbative expansion with respect to 
the coup ling constant. Examples for the first order of this expansion are considered in a number of papers, see for 
example (la. H3|. 

The T-matrix representation for the electromagnetic field has the same general structure as that for the scalar field. 
The main difference is in the presence of the polarizations. The expansion basis for the electromagnetic field has two 
components, 

m Mm (r) = u M ; m (r) = L— ^= u k .i m (r) , (16) 

v L 

nfe.im(r) = u fc , 2 im(r) = V x L = u k ^ m {r) , 

V — A vL 

where L is the orbital momentum operator. Under a translation these functions mix and in place of {?]) the translation 
formula is now 

m Mm (r + ae z ) (17) 
= E ( B im,i'm> (a) mk,i'm' (r) + Ci m ,i' m > (a) n k< i> m > (r)) , 

I'm' 

n fe ,i m (r + ae z ) = 

= E ( c im,i'm'(a) m fei ;/ m /(r) + B lm ^ m >(a) n fe j/ m /(r)) . 

I'm' 

The Greens function (now in fact the Greens dyadic) can be expressed in the basis (|16l) similar to and reads 



2 f°° dk k sr^ 
G f (r,r') = - / 2 + fc2 2^ Ufc,^m(r)i4, s ; m (r'), 

shn 



(18) 



s taking values s = 1,2 and the s-wave is excluded, I > 1. The translation coefficients Bi mi it m i(a) and C; m! ;' m '(a) are 
known in the electromagnetic theory, as pointed out in [l5[ a particularly useful representation can be found in [T^ |. 
In the given geometry, the coefficient -B; m ,i'm' can be obtained by the substitution 



< - «; (19) 
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in A lm , Vm , with 

A r = \ [l"{l" + l)-Kl + l)-V{l' + l)\ 

y/l(l + l)l'(l' + l) 

Cl m ,l'm' — A^' Ai m j< m i (21) 

A U , = 2m ^ d + R) (22) 

y/l(l+l)l'(l'+l)' 

where the Ai m ^ m i are the same as in the scalar case, Eq.©. The coefficients A l u , and Aw result from the orbital 
momentum operators, for instance from the normalization factors in (|16p . 

The translation from one coordinate system to the other mixes the polarizations of the electromagnetic field. 
Therefore the energy has now the representation 



The other one is given by 



with 



E = - [ ^Trln(l-N). 
2 7_oo 2^ 1 ' 



(23) 



and the trace is 



Tr= E tr < ( 24 ) 

m=0 l=max(l,\m\) 



where tr denotes the trace over N which is now a (2x2)-matrix in the polarizations, 

l+V 



Y, (25) 
^ l»=\l-l>\ 

( ai;;, h,v \ ( df B (£R) o 

K [A ltl , A? tl ,)\ -d[ M ((R) 

The functions df E and df M describe the scattering of the corresponding polarizations of the electromagnetic field on 
a conducting sphere and are similar to that of the scalar field. For the TE mode it is literally the same, 

<neB) = Z +1/2 Tl . (^) 

The minus sign in front of df in (|25[) results from the the spin of the electromagnetic field under reflection on the 
plane. 

Representation (|23|) of the vacuum energy of the electromagnetic field was derived in different notations in [1, [t| . The 
coincidence with these formulas can be checked by comparing some first orders of the expansion for large separation. 



and for the TM-mode it is 



III. THE ASYMPTOTIC EXPANSION FOR THE SCALAR FIELD 

In this section we consider the asymptotic expansion for the scalar field at small separation. We follow [B[ and 
repeat the main steps of the derivation since these appear essentially in the same form in the next section for the 
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electromagnetic case. At once we add a discussion on all combinations of the boundary conditions on the plane and 
on the sphere. 

The vacuum energy is given by Eq. (fl5| . First of all we make the substitution £ — > £/R to get rid of the dimensional 
variables. Then, as already mentioned, for decreasing separation the convergence of the integral and the sums in (|15p 
slows down and the main contribution comes from higher and higher frequencies and orbital momenta. We know, 
by hindsight, the region delivering the dominating contributions. Since at small separation all summation indices 
involved take high values we substitute all sums by corresponding integrations. In this way we drop exponentially small 
contributions which is allowed aiming for an asymptotic expansion. In these integrations we make the substitutions 




e 

(28) 

where we divided the orbital momenta by means of li = I + li (i = 1, . . . , s) into the index I of the main diagonal 
and the off-diagonal indices li. The variable r has the meaning of the cosine of the polar angle in the £, ^-plane. This 
substitution describes the region where the main contributions come from. 
In the new variables the expression for the energy reads 

E = - AV— dtte- 2t ^ 

(29) 




where 



"o^-W^ + ^HV^ (30) 



is the lower boundary in the -integrations. It follows with (|28[) from I > \m\. In 



Z -11^ — ^+^^' ^ 

collects the information from the scattering process together with the prefactors which follow from the substitution 
(|28p. In (|31[) we use the formal definitions Iq = l s +i = and we defined 

s + 1 

r7 as = 2t(s + l) + r/ 1 +^ 2 (32) 

T 

with 771 = J2i=o( n i _ n i+i) 2 - 

Next we expand Z for small e. It turns out that it is possible to do this expansion straightforwardly by expanding 
all quantities entering. These are the Bessel functions and the Clebsch-Gordan coefficients in the iVj j/, Eq. fTQ]) . Here 
we can follow the corresponding expansion in While for the Bessel functions the known uniform asymptotic 
expansion can be used, for the Clebsch-Gordan coefficients the corresponding asymptotic expansion was first derived 
in It rests on an integral representation of these coefficients. In this way, one arrives at an asymptotic expansion, 
Ni,v with 

OO 

= / £T e -2t-(n-n>f f *^ -^M^W" ( 33 ) 

—00 

!/=0 ■ V ' / 
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The function f(n, u, tr) and the corresponding function in the next order can be found in The integration over n 
results from the mentioned integral representation and the sum over v from the summation over I" in (|10[) with the 
substitution I" = I + I' — 2v. By the symmetry properties of the 3j-symbols in (jX3j) it follows that v takes integer 
values only. The upper limit of the summation over v is 

u m = \{l + l' (34) 

The exponential factor follows from the function n{z) in the asymptotic expansion of the Bessel functions. 
By means of the substitution ([28]) . the upper limit of the v summation depends on e, 

v m = \ — \n - n \. (35) 

eye 



In the considered case of a scalar field it is possible to put e = here. After that the sum over v in (1331) can be carried 
out. The integral over n is Gaussian and can be carried out too. We get 

W# (36) 

= ^ ^-^'^ + <%? >Te + a% e + ...). 

The functions a^/ 2 ^(n,n') and a^(n,n') are given in Eq.(A.22) in Q. 

Eq. ([3^|) must be inserted into Eq. ([3T|) . The prefactors and the exponentials just cancel and the remaining depen- 
dence on e is contained in the bracket. This fact justifies the substitution When inserting N^, into (|3"Tj) we get 
the asymptotic expansion Z as . After a re-expansion it takes the form 

Z" = l + T l ^ / ± 1 Vi (37) 



i=i 



+a D e + 



with 



«°= E ^L^t+E^nr (38) 

0<i<j<s i=l 

We mention that we used in the re-expansion the general formula 

s s 

Y[(l + Xi) =l + ^2xi+ x i x j + ... 1 (39) 

i=0 i=0 0<i<j<s 

which holds in the sense of an expansion for small Xi. We note that the product turned into sums. We will use this 
formula below several times without further notice. 

Inserting this expression for iJ as into ([29")) we get the asymptotic expansion of the energy. Here we have still an 
e-dependence in no, (|43[) . However, we are allowed to put this e — 0, i.e., to take no = — oo, since in doing so all 
integrations remain finite. After that the integrations can be carried out. These are either over simple exponentials 
or are Gaussian. We mention the most complicated one, which is that over the nj. It was calculated in Q, Eq.(66), 




The result for the energy is 



— r- (40) 

Vs + 1 



Carrying out the summation over s, in the leading order just the PFA emerges. In the order y/e the integrations gave 
zero for symmetry reasons and in order e we get the first correction beyond PFA. 
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The energy, Eq. (|4Tj) . is for Dirichlet boundary conditions on both surfaces. The case of Neumann boundary 
conditions on the sphere can be handled in complete analogy. There are only two differences. The first one is an 
additional minus sign resulting from the derivative of the Bessel function K v in the denominator in (fT2f . It appears 
in each factor Niji, hence it gives a sign factor (— l) s+1 to the sum over s. At this place we restore the notation of 
Eq. (fT4]) for the different combination of the boundary conditions and accounting for all signs we get 

xn 1 R ^ / i \ 

XN l ~ 

16tt d 2 ^ (s + l) 4 

s— v ' 

*(i+(H<.+i>" 

The second difference for Neumann boundary conditions on the sphere is an additional contribution containing the 
factor of (s + l) 2 in the last line. It results from the change in the Debye polynomials due to the derivatives in (JT3J) 
in the asymptotic expansion of the Bessel functions and results in a different function a N in place of (|38p . As shown 
in this factor is the same even if generalizing to Robin boundary conditions. For instance, it is the same for the 
TM mode of the electromagnetic field in the next section. 

In Eq. (|42| the summations result in Riemann zeta functions; we need 



oo 

y L 

^ (s + l 



oo 1 n 4 °° (-1) S 7 

+ =C(4)= 90' (s + I) 4 = 8 C(4) ' 

s—0 s— 

In each case the leading order in (142(1 delivers the PFA. The relative corrections are 

£ DD 1 

! + -£+..., 



EppA 3 

£ NN /I 10 

1 



E PFA \3 it 2 

(\ 5 



1 



r'2 



EpFA \3 7T 

= 1 + - £ + ... . (43) 

We note that the corrections for Dirichlet boundary conditions on the sphere, but Dirichlet or Neumann boundary 
conditions on the plane (first and last line) are the same. This follows simply from the structure of the signs in (|42[) . 



IV. THE ASYMPTOTIC EXPANSION FOR THE ELECTROMAGNETIC FIELD 

The asymptotic expansion for the electromagnetic field start with the same steps as in the scalar case, i.e., with the 
substitution of the orbital momentum sums by integrals and the substitution J28|) . The next step is the asymptotic 

expansion of the functions N^f ' entering the matrix N; ;' in Eq. ((25|) . In the functions dJ E , Eq. ([26|) . and d™, Eq. ([27|) . 

we use again the uniform asymptotic expansion of the Bessel functions. In the factors A l w and A l w , Eq. ([2H]) . we use 
directly the substitution (|28[) and write them in the form 

A%, =l + e\ n , n > (44) 

= 1+ ( _ A _ 2™ n{n)+n{n') + e 2 v(2v - 1) 
\^(n)j(ri) J y/t j(n)j(n') t 7(^)7(71') 
= l + A + B + C, 
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7W7F) 



with the notations 



= tV< + 2nVe, 7(71) = \/^(n)(/x(n) + e/Vt) . (46) 

Here /x(n) follows from Z and 7(71) from / + 1. We divided A' j, into three parts which will be treated separately in the 
subsections below. All these will deliver a contribution of order e (including elne and e(ln£r) 2 ) although this cannot 
be seen directly from Eq. (|44|) . The same holds for A^s Eq. (|45|) , after taking the trace over the polarizations. 

The energy is given by Eq. P5|) and its asymptotic expansion by Eq. (pj"9")) with another Z which is, of course, still 
similar to (f3Tj) . It is expressed in terms of the matrix Ni t i>, Eq. ([25|) . by 

For a matrix N; p we get the asymptotic expansion 







r-Vi ( 


(0 







,(|)(TE) r \ / (1)(TE) f 



Here the functions a n ,n' are the same as in the scalar cases with the corresponding boundary conditions. The factors 
in front of the figure bracket appear in the same way as in the scalar case from the Bessel functions. We remind that 
in the given order of the asymptotic expansion the difference between scalar Neumann boundary conditions and those 
of the TM mode does not show up. 

In deriving Eq. (|48|> we have to pay attention to the t'-dependence of \ n ,n' ■ Therefore we defined \ n ,n' by 



00 v 



Vn'. ( 4 9) 



This definition is taken in a way that A„, n ' — > 1 in (|49p gives A n>n i = 1. 

Next we have to insert (|48l) into (I47[) . Making a re-expansion in e and taking care of the matrix multiplication we 
get 





E 

0<i<j<s 



(50) 



In this formula we dropped off-diagonal contributions and those proportional to \fe which will disappear later when 
carrying out the remaining integrations as already mentioned in the preceding section. In this formula, the factors 
resulting from the T-matrix collected in the same way as in the scalar case into the functions aP and a N . The 
contributions from the vector structures (|44[) and ([43]) do not depend on the polarization (see the structure of the 
translation formulas, Eq. (|17[) ). Therefore these enter proportional to a unit matrix. Now we take the remaining 
trace and come to 

Z as = 2+ ( a D +a N + 2^A n! ,„ !+1 +2 £ A ni) „ i+1 X nj ,n j+1 ] e + • • • (51) 

V i=0 0<i<j<s I 

This expression for Z as must be inserted for Z into the energy, Eq. (j29f . At this place we can already read off some of 
the features. First of all, the factor of 2 accounts for the polarizations of the electromagnetic field. Since the remaining 
integrations are the same as for the scalar field we immediately get the expected result that in PEA the electromagnetic 
field gives twice the contribution of a scalar field. This is the leading order and it is of course independent of the 
boundary conditions. 

In order e, the first two factors, a D and a N , give the same contributions as in the scalar cases. So we have, for 
the electromagnetic field, in the first correction beyond PFA a contribution with Dirichlet boundary conditions on 
both, the sphere and the plane, and another one with Neumann boundary conditions instead. The corresponding 
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contributions to the relative corrections beyond PFA to the energy are given by the first two lines in Eq. (|43|) and 
must be divided by 2. 

We mention the role of the minus sign in the latter case. For the scalar field there were two of them, one following 
from the derivative of the Bessel function in the denominator in (|12[) . This sign has a corresponding one in the 
electromagnetic case in the denominator of (I27p . The other minus sign followed in the scalar case from the sign in 
the logarithm in Eq. (fT5|) . In the electromagnetic case the corresponding one is the minus sign in front of dJ M in (|25|) . 

The remaining terms in Eq. (|5ip just represent the additional contributions which come in from the vector character 
of the electromagnetic field, i.e., from its spin. There are contributions diagonal in the polarizations, resulting from 
A\ Eq. ([20|) . and off-diagonal ones resulting from A;,;/, Eq. (j22)) . It should be mentioned that in the considered first 
order in e all these contributions enter additivly. A mixing of these will happen in higher orders only. 

In the following subsections we consider separately the contributions resulting from the three parts, A, B and C 
of K\ t ,, Eq. (|4"4"|) and that of A;.;*, Eq. (|4"5)l . Starting from part B we will meet expressions where it is not possible to 
make a direct expansion for small e. For illustration we consider a simple example. Consider the integral 

1 

f{e) = [ dr g(r) - „ (52) 

y ' J r + ae + be 2 K ' 

o 

for e — > 0. In case the function g(r) has a zero, g(0) = 0, we can put e = directly under the sign of the integral. 
In opposite, if g(0) ^= holds, we cannot do that since the r-integration would diverge. The only we can do is to put 
e = where it goes with the coefficient 6, 

l 

/(e) = [dr^- + ..., (53) 
J t + ae 
o 

where the dots denote contributions of higher order in s. The remaining integral must be treated in some other way. 
For example, we can integrate by parts and expand after that, 

f(e) = ln(l+e) 5 (l)-lne$(0)- / rfrln(r + e)g' (r) + . . . , 

Jo 

= - In e.g(O) - dr lnr g' {t) + ... . (54) 
Jo 

Below such and similar situations will appear repeatedly. 



A. Part A in the A-contribution 



Part A is given by 



A=^l-l (55) 
7(n)7(n'J 



with n(n) and 7(n) defined in Eq. (|46|) . The contribution from this part to the energy is quite easy to calculate since 
it is possible to expand A directly in powers of s, 

A = ~e + ..., (56) 

IT 

without causing any divergences. We define its contribution to \ n>n i by 

^A;n,n' = ~ ~- (57) 
tT 

This is a quite simple formula, for instance the dependence on n and n' dropped out. Further, since it does not 
depend on u, from (|4*5|) we have \A;n,n> — ^A;n,n' an d its contribution to Z as , (f5"T) . is 
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Here the dots stand for all other contributions. We denote the corresponding part of the energy by AEa and we come 
with dm to 

R —2 f°° f 1 dr \fr f°° dfi 



R ^ -2 f , f drJr f°° dn 



X 



\i =1 



Y[ I ^ ] ( 1 _ e £±1 + . . . ) e -2t( S +l)-M 2 ( S +l)/r-r ;i _ (5Q) 



We mention that we have taken the lower integration limit in the n 3 ^-integrations to — oo which is possible since the 
remaining integrations do converge. The integrations in this expression can be carried out easily and the result is 

E PFA + AE A = — l__ £ — + ... 1 . (60) 



The correction resulting from part A is 



B. Part B in the A-contribution 

Part B is given by 

B = -2ve n(n) + /x(n') . 
y/i 7(71)7(71') 

Regrettably, its contribution to the energy cannot be calculated so easy as before. First of all we have an additional 
dependence on v. In addition, it is impossible to make a simple expansion in e. This would produce a singularity in 
the T-integration. 

We define \B;n,n' = B/s which must be inserted into (|49p . The summation over v is quite simple, 

Here we have taken u m — oo since this does not cause singularities. Next we need to carry out the integration over 77 
in It is Gaussian, 

00 

/ ^= * exp (lT7 7/2 ~ 7/2 + 2lv ^ + /i2 ) = Mt ' m) V l^T cxp ( _ ^) ' (64) 



l-T /, „1 + T 2 



with 

/ lB ( r , M ) = 2_M-2^— ^1. (65) 
Using Ea.([M| in |gS} we get 

As ; „, „' = fts (t, /«) As ; n,n' ■ (66) 
This must be inserted into i? as , (fSTj) . and further into the energy. With (f2"5]) it is 

where we defined 

s= r ^ e-^^ hfl (r^) f n r^E M ? ) w (T ^ l) e " m - (68) 
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The integration over \i is Gaussian and we get 



t4jH) + fJ,(n i+1 ) 

H 7(^)7(^+1) 



(69) 



with 



1 /1-T / _ 1+T 

4V 1 + t 1 s + 1 



(70) 



Now we have still e in a number of places in (|69|) . Still we cannot simply expand for small e. If we would do so, 
because of 



1 



(71) 



where we used (|46|) . the r-integration would become logarithmically divergent. The only places where we can put 
e = directly are those where it goes with n in //(n) and in 7(71), 



1 



(72) 



Also we can take no = —00. In doing so we do not produce divergences in the integrations. After that we are left 
with a simpler r-integration. Here we integrate by parts, 



Jo 



(It 



Ib{t) 



dr In 



VrVt + ^Ty/i + e/y/i g 



Ot 



fair). 



(73) 



The surface term is zero. In the new r-integral it is possible to expand the logarithm for small e. We insert the result 
into B, Eq.flBSJ), and get 



B 



-hn£ + hnt + In 2) (f B (l) - f B (0)) + 1 



dr lni 



d_ 

dr 



f B (r) 



(74) 



Since after (|72p there is no more any n^-dependence we used formula (I40[) and accounted also for the sum over i. The 
remaining integration over r can be carried out easily, 



1 d it l + sln2 

dr lnr— f B (r) 



di 



8 4(s + l) 



(75) 



We mention that Eq. (|74[) is the first place where a logarithm in e appears. Its origin is clearly seen from Eq. (|7ip . 

Next we have to insert B into the energy (I67p . The remaining i-integration is now a bit more complicated since it 
involves a hit. From that a logarithm ln(l + s) appears. However, simple calculations yield 



R 
And 2 



(C(3) - C(2)) (7 - 2 In 2) + C'(2) - C(3) - \ C(2) + C(3) + i- (C(3) - C(2)) lne 



(76) 



This is the contribution from part B in (I43f> to the corrections beyond PFA. It involves a logarithm in e and it has an 
analytic expression. 



Part C in the A-contribution 



Part C is given by 



C 



£ 2 u{2v- 1) 
t j(n)j(n') 



(77) 
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The calculation of its contribution to the energy requires most effort. First of all we observe a quadratic dependence 
on v. We have to insert (|TT|) into (f3"3"|) and define the corresponding contribution to Nffi by 



N, 



C;l,V 



ST 



2vrf(l + r) 



-2t-(n-n') 2 f drj r/ 2 +2ir;V2^+^ 2 



(78) 



v=0 



A VI 



C. 



Here it is impossible to take the upper limit v m of the summation over v to infinity. Doing so would produce a factor 
t -2 and making the r-integration diverge. Therefore we must account for a finite u m , Eq. (|34p . With the substitution 
it becomes 



tT 



(79) 



Technically we account for it by inserting a step function into the sum and formally summing up to infinity as before. 
For the step function we take the integral representation 



®{v m -v) = 



dp exp (i p (±L-y)) 
2m p — iO 



and define 



N, 



C:l,V 



dp exp(ipf) - 



N, 



C:l,V 



I _ ca 2m p — iO 
The summation over v now involves the additional factor a" with 

a ee e~ ip . 

This sum and also the integration over r\ can be done generalizing l|63p and (164p . The result is 



N, 



C;l,V 



t 7(71)7(71') 



-77„-^ a /(TA) 



4?rf 



(80) 



(81) 



(82) 



(83) 



with 



hc(r,fJ,) = 6 



1 - T 



4r 

1 — T 



A 



1-4 A/i 2 + - A /1' 

r 3 \ t 



4r 



a A 



1 - 2- 



1 + t 



A/x 2 



Here we defined (for use only in this subsection) 

2r 



A : 



l + r-a(l -t)' 



A 



2-(l-r)(l-a) 



(84) 



(85) 



From (f8"Tj) we can now read off the contribution from part C into A„ jra < and insert that into Z as , Eq. (f50")) and further 
into the energy (|29[) . The corresponding contribution is 



AEr 



R 

Ami 2 



00 9 Jl poo 



(86) 



where we defined 
- 1 



lo vT 



-/x 2 (s+l)/r 



n 



E 



t / ~^ l{nih(n i+ i) J-oc 2m p-iO 



dp exp (ip^-) 



hc(r,n). (87) 
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Here, again, the integration over /i is Gaussian and we come to 



C = 



I f 1 dr t I W f°° drij \ e Vl 



Vs + l Jo Vl 



n 3 



i?(r,i) 



with 



and 



fl(r,*)= / tt£? 



dp exp(ip^) / -^1/2 



27ri p — iO 



AA 



6(aA/ c (r)) 2 + aA/ c (r; 



(88) 



(89) 



(90) 



In Eq. ([88|) . the main contribution comes from r ~ 0. This is because the function fc(j) diverges like fc( T ) ~ r 
for r — * 0. The integration over r is nevertheless finite for e ^ because of the exponential involving p. Next we 
want to carry out the integration over p in (|89p . For this we move the integration contour upwards in the complex 
p-plane. We have a pole at p = 0. Its contribution is easy to calculate and it gives just the result we would obtain 
with putting v m = oo at the very beginning. Further there are poles in po = i In y±^ + 27m {n integer) from A and 

cuts starting from p c = i In + ir + 27m resulting from AA. 

For small e, non- vanishing contributions come only from the poles in p — and in p = po with n — 0. For these we 
can expand po — 2ir + . . . and AA = 1 + . . . and R takes the form 



dp exp (ip^-) 
2ni p — iO 



6 a 



The pole in po is second order and we get 



R(T,t) = 6f c (T) 2 



1-1 



2ir 2 



P-Po 



-2tr 2 /e 



Mr) 



-Po 



P~P0 



fc(r) 



+ fc(r) 



1-e 



-2tr 2 /e 



(91) 



(92) 



We mention that the first term in both square brackets result from a ^-summation with v m = oo. The terms with the 
exponentials result from taking a finite v m . As a result, while e ^ 0, the function J2(r, t) is finite for r — > instead of 
diverging like r" 2 as the function /c(t) does. We emphasize that this is the justification for accounting for a finite 
v m in part C. In the other parts, because we did not have a small-r behavior like in this one, we could take v m = oo 
without hitting a divergence. 

The expression for R, Eq. ([9^)l . must be inserted into C, Eg. ([55]). To proceed with the expansion for small e we 
observe that we cannot put e — in R since this would return us to a situation where the r-integration diverges at 
t — > 0. Instead we make in C the substitution t — * T\fe. 



C=~ 



1 1 




R(Ty/e,t) 



with 



and 



7(n) = ^7 



rVt + 2n) \rVi+2n 



n Q = --rVi. 



(93) 



(94) 



(95) 



In this way, the r-integration produces a factor 1/e which makes the contribution from part C, which initially went 
with a factor e 2 , a contribution first order in e. We mention that it resulted from the factors j(n) in the denominator. 

There is still a dependence on e in C, Eq. (|53")) . It is twofold. First is that which goes with r. Here we can put 
e = 0. The same we can do in the upper integration limit. We note 



g(2tT 2 ) + 0(l) 



(96) 
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with 



for £ — ► and for C we get 



l(x) = \{l-{l+x)e- x ) (97) 



~ 3 s 2 f°° dr . ». / A [°° dnj \ A 



In order to simplify the representation of C we make in (1981) the substitution r — > t \Jt. After that it takes the form 



C-./ ^)MI/ ^l^^-j. (99) 



where we introduced the notations 



3 s 2 r 



" 2e( s + 1)5/2' "o 2 , 

which will be used in the remaining part of this subsection. 

In this way we are left with the dependence on e in 7. Here we can not put directly e = 0. This would produce 
a divergence in the integrations at r = — 2n and an imaginary part would appear which is clearly not present in the 
energy. The way out is a partial integration in the ni and in the n^+i integrations. But before we can do that we 
have to pay attention to the contributions from i — and from i = s in the sum over i. From the formal setting in 
Eq.(|3ip we have to put jiq = n s = in Eq. (l9"9"l) and the corresponding 7 do not depend on any n. We denote the 
contributions from i = and i = s (both give for symmetry reasons the same contribution) by Co and the remaining 
one, i.e., that for i = 1, . . . , s — 1, by C\. 

First we consider Cq. Since the function g(x) ~ x 2 for x — > we can take 7(0) = r + . . . and are left with 



7(n) = \]{t + 2n){r + 2n + y/IJi) , (100) 



Jo T \j=V ft ° / 7( ni ' 



r 'n = ^ / -.</<-^) Ml/ ^ 1 77T^- (101) 



Now we integrate by parts according to 



e 



-m ( / £ \i/4 



dn — — = - In I 2 

r/2 7W 



(?) 



-r/2 



9 



with 



Here we can expand for small e, 



dn In ( 2 ( Vr + 2n + + 2n + i/s/t) ) ) — (102) 



m = m\n=- T /2- (103) 



p-'/i / /P\l/4\ f 00 f) 

dn— = -In 2 - } e~ m - dn In yr + 2n) — e^ 1 + ... . (104) 

r/2 7W V Vi/ / ./-r/2 ™ 



This must be inserted into Eq. (|101[) . After that we make there the substitution nj — > nj t/2. This allows to change 
the orders of integrations, 

C = ^ filn§+]n2) i?(s) + 4ctQ(s) , (105) 
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and we introduced the notations 



R(s) 
Q(s) 



nf£i»(? 



i-nz)g 2 — 



(106) 



with 



3i (^) 
92 0) 



— 5 



T 
72 



<*">(§)' 



<?(2r 2 ) In (A^^l + 2m) 



s+l 



(107) 



The integrations over r can be carried out explicitly, however the formulas are too voluminous as to be displayed 
here. The remaining integrations over the nj can be performed only numerically. Even that is not an easy task since 
the integrals are s-dimensional. We could proceed only till s — 5 for R(s) and s = 3 for Q(s). However, since the sum 
over s is quite fast converging this gives at least the order of magnitude correctly. 

Expression (|105[) must be inserted into the energy (|86[) . Introducing the corresponding notation it is 



AE c 



R 

And 2 



00 o c-2 /- 00 

Y^L_ / dte- 2t ^C . 



(108) 



Here, the t-integration can be carried out easily and summing over 5, as far as data are available, the result is 

AE 6q = -ip£ (0.0020 + 0.00017 hie) . 



(109) 



Now we have to consider Cx, i.e., the contributions from i — l,...,s — 1 to (|9"9")l . We make the substitution 



rij t/2 such that it takes the form 



C x =o\ $ 9 (20(- 



a -r;iT 2 /4 



71 ) 1~[ 7("t)7K+l) 



(110) 



with j(n) = \/ (1 + n)(r(l + n) + We/t. For the integration by parts we adopt the following scheme, 



e -»)iT 2 /4 

dn — — r-r — 

i l{n) 



- In - - In (2y/r) 
4 t v ' 



-??2T 2 /4 



oo V?Vl + n + \/t(1 + n) + VeA 9 

drain V --" 1T 



-171 t/4 



-1 



/r • 



<9n 



(111) 



Here we can take e — > and get 



-r;iT 2 /4 



1 l(n) 



1 f'°° 

4= / dn (Li (») + L 2 (n)) e -"i- 2 /4 + 



where we defined 



1, e 



L 1 (n)= ( 7 ln 7 -ln(2V7) )5(n+l), L 2 (n) = In VTT^ — . 



4 t 

We have to apply these formulas in (|110[) to m and to nj+i, 
Ci = 



/*oo J I ^ poo J \ s — 1 

a / 4 3 (2r 2 ) Q' A / -4 E + L2 ^)) (iiK+i) + ^(ni+x)) e""^ 



/4 



(112) 



(113) 



(114) 
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Multiplying out the two brackets, 

(Li(rii) + L 2 (rii)) (Li(n i+ i) + L 2 (n i+1 )) = Lx{ni)Li(n i+ i) + 2L 1 (n i )L 2 (n i+1 ) + L 2 {ni)L 2 (n i+1 ) (115) 
(we used the symmetry under i — > s — 1 — i in (|114j) ). we split C\ into three parts, 

C 1 = C 1A + C 1B + C w (116) 

and consider them separately. 

We start with C\a and interchanging the orders of integration it is 



Gl4 = (77T ' 



n / dn 3 W 5 3 (2^)(l) S (iln£-ln(2^)) (117) 
3=1 J 1 / i=i J ° V / 



with 

% = »7iK= ni+1 =-i- (US) 

This is the place where the logarithm of e appears squared. Again, the integration over r can be carried out explicitly 
delivering lengthy formulas. Expression \\YJ\ must be inserted into the energy (|86[) . Introducing the corresponding 
notation it is 

2 00 F 2 

A ^ = -^E7Ti (T7r ~ s/2p(s) ' (119) 

where 



/» CO 

/ ^e- 2 *( s+1 )(7 1A = ^- s / 2 P( S ) 
Jo 



(120) 



collects the integrations over the 7ij and over t. Again, the n-integrations must be done numerically, we were able to 
go up to s = 7. As a result we get 

AEcia = ^£(-8.810- 7 -2.410~ 7 lne-3.61Cr 7 (ln£) 2 ) . (121) 



-4a.-/ 2 f f[ jTd^j £ jf° ^^r 2 ) (I)** Q In f - In (2^)) (2n J+1 



Next we have to consider C\b- In parallel to (|117p it is 

(122) 
with 

V2 =m\m=-v (I 23 ) 

Again, the integration over r can be carried out explicitly and the integration over the rij only numerically. Here we 
could go until s = 5. The result is 

R 

AEcib = -TT £ (-9.4 10~ 6 + 0.000019 Ine) . (124) 
dr 

Finally we come to C\c- In parallel to (|117p it is 

JiJ l dn AllJ o ^^(2r 2 )(I) S ln(l + n 4 )ln(l + n l+1 ) 

x ^1 + y(2n 4 - - n i+ i)(2n J+ i - n< - n l+2 )^ e^ 11 " 2 / 4 . (125) 

As before, the r-integration can be done explicitly and the n-integrations numerically. Here we went up to s = 5. 
The result is 

AEcic = 5 £ (-0-000076) . (126) 
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D. The A-contribution 



We start from the formula (|45|) and insert X n , n ' into -Z as , Eq. (|5Tj) . From this we define with ([29]) the corresponding 
contribution to the energy, 



oo poo 

J dtte-^+^L (127) 



s=0 



with 



£ = /' ^ p *L c -M"(-+i)/r (fl r*±\ y 4rt(l-^)(l- £ )V ^ ( | . )s . 

'o \/l - t 2 i_oo 0? V=V™o J a< ~^ J<s T(n.)7(%i)7(nj)7(«j+i) 



Integration over /z delivers 



i- 1 r d jn f]r dn A t ^a-- 2 )(i-^ c - m (129) 

(s + lWo >/T^ \JLV»o V^ y | 0< ^ <s 7(n i )7(«m)7K-)7K- +1 ) ' 1 j 

Now we can put e = where it goes with rij and in uq, again keeping all integrations and summations finite. With 
this, we note 7(71) = Vrt\/ t + e/t + . . . . As a consequence, the rij-integrations become simple and are reduced to 
formula (|40p . Also the dependence on i disappears and the sum is X)o<i<j<s = s ( s + -0/2- I n L only one integration 
remains, 

1 ' r -\/l — r 2 s ( t 



L = r / dr — T = h - . (130) 

t(s + l)J ( r+f ) 2 t(s+l) \ej { ' 

This expression must be inserted into the energy (|127[) . Taking into account h(t) = ln(2i) — 2 + . . . for t — ► 00 we can 
carry out the ^integration, 

^ = 53! -t 77^ (-5 " 1 " I " 5 1,1(8 + ") + ' ' ' ' (131) 



a + 

Finally we carry out the summation and come to 

i? 



(C(2) - C(3)) (2 + 7 ) + C'(3) - C'(2) + ^ (C(2) - C(3)) hxt 



V. CONCLUSIONS 



(132) 



In the forgoing section we calculated separately the different parts contributing to the correction beyond PFA in 
the electromagnetic case. These are the 'scalar' ones (from the first two lines in l|43p). 

AE TE + AE TM 15 , 
^ — = ---5«-0.173, 133 

that from parts the A, ([6T|) . and B, (|76|) . in the Lambda-contribution and that from the A contribution, (|132p . 



AE A + AE B + AE K _ 180 
eE PFA 7T 4 



[(1 + 2 In 2)C(3) - 2(1 + In 2)C(2)] « -4.99 , (134) 



and that from part C in the Lambda-contribution, Eqs. ([TD9]) . (fl"2"Tj) . (fl"24|) . ([T26]) . 

+ A^cia + A£ C ib + A£:cic 



eE'PFA 



0.045 - 0.00441ne + 8.5 MT^lne) 2 . (135) 



In (|134|) cancellations happened, for instance the logarithmic contributions compensated each other. This contribution 
could be calculated analytically, like the 'scalar' one, Eq. (|133[) . The remaining contribution (|135p could be calculated 
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only numerically. It contains the remaining logarithm and it is numerically small. This smallness justifies the use 
of the quite small precision reached in subsection C. Putting all together we come to the relative correction for the 
electromagnetic case, 

AE EU 

— = -5.2 - 0.0044 lne + 8.5 10" 6 (lne) 2 , (136) 

from which representation ([3]) follows. The main contribution is a constant, the logarithmic contributions are numer- 
ically very small and do not play a role at any reasonable separation. 

The result (|136[) is quite unexpected since it it quite large. It must be mentioned that it is not in agreement with 
the numerical calculations performed in [1, Q, where ~ 1.4 for the constant contribution was obtained. In these 
calculations the fits were made without accounting for possible logarithmic contributions. However, in view of the 
smallness of the coefficients in (I136[) this should be acceptable. 

The result (|136[) also does not support the experimental results found in ((J. However, the experiments were done 
with real metals. It should be a subject of future work to account for that in an analytical calculation. 

In general, it must be underlined that Q136P is the second term in an asymptotic expansion. Therefore it is hard 
to predict how small e must be to get a good approximation in this way. In principle, it cannot be excluded that 
(|136[) gives a good approximation only for e smaller than that which are interesting for the experiments and which 
are accessible by the mentioned numerical approaches. 

There is still another problem with the first correction beyond PFA for Neumann boundary conditions. Already 
in the easier case of a cylinder in front of a plane the numerical calculations reported in [jj showed good agreement 
with the analytical ones in Q only for Dirichlet boundary conditions, but not for Neumann ones. For the latter, 
only with a fit including a logarithmic term in the second order, e 2 lne, agreement was found. For a sphere in front 
of a plane for a scalar field, the numerical results reported in Q are in agreement with the analytical ones, Eq . (|4"3")l . 
only for Dirichlet conditions on both, the sphere and the plane (DD). In the other three combinations of boundary 
conditions the numbers are quite different. This is quite unexpected for the case of Neumann conditions on the plane 
but Dirichlet conditions on the sphere (ND), since the analytical results for (DD) and (ND) are the same. 

It should be mentioned that the appearance of logarithms in the expansion is probably a rather common feature. 
This can be seen from the structure of the corrections, see, for example, Eqn.(B13) in Q. The expansion parameter 
e is always accompanied by a factor 1/t, producing in higher orders a singularity at t — > 0. It is only in the order e 
considered in jij as well as in Q that these did not show up. 

In view of the agreement of the numerical results with the analytical results in the (DD) case, the disagreement 
in the other three cases where we have at least on one surface Neumann conditions can be viewed as a hint that the 
numerical approach is more difficult once Neumann conditions are involved. It could happen that an agreement can 
be reached for smaller e only. We would like to point out that also in the analytical approach the calculations with 
Neumann conditions are a bit more delicate. The point is in the convergence of the sum over s in (|29p. While for 
Dirichlet conditions the decrease is ~ (s + 1)~ 4 (first line in I02J)), it is only ~ (s + 1)~ 2 for Neumann conditions. In 
the electromagnetic case it is even weaker, ~ (s + 1)~ 2 ln(s + 1), for example in (|13ip . 

As a consequence of the mentioned disagreement it would be interesting to improve the numerical approach. The 
main obstacle is that very large orbital momenta must be accounted for. A way out could consist of three steps. First, 
one may expand the logarithm as in Eq.fll5p. As we know from the analytical approach this sum is converging. To get 
a satisfactory precision, to take a few terms should be sufficient. In the second step one would make the substitution 
(|28p also in the numerical approach. This allows to capture the main contribution. Finally, as third step, one would 
need to adopt some approach of coarsening to the orbital momentum summations or their substitution by integrals. 
In any case, further work is necessary in this direction. 
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